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DEFINITION OF SYMBOLS 



FIGURE 1. CONFIGURATION OF A VIBRATING CANTILEVER BEAM 
Symbol Definition 

x axial coordinate of beam (in. )* 

L beam length (in.) 

z = x/L dimensionless axial coordinate of beam 

y lateral deflection of beam (in. ) 

Y Q initial displacement of free end of beam in free vibration (in. ) (See top of 

next page for definition in forced vibration) . 

* These are the most common set of units used in the U.S. A. ; however, any consistent 
set may be used as the formulas are presented in dimensionless form. 
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Symbol 

Y 0 
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co 
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DEFINITION OF SYMBOLS (Concluded) 

Definition 

amplitude of lateral displacement of "fixed" end of beam in forced vibration 
( in. ) ( See previous page for definition in free vibration) . 

weight of beam (lb) 

moment of inertia of beam cross-section (in. 4 ) 

modulus of elasticity of the material of the beam (lb/ in. 2 ) 

acceleration of gravity, (in. / sec 2 ) (taken equal to 386. 4 in. / sec 2 in sample 
examples) 


Elg 


j, beam parameter in free vibration (sec -1 ) 


WL 


circular frequency of excitomotor (rad/ sec) 


WL 3 O J 2 
Elg 


% 




beam parameter in forced vibration (dimensionless) 


coefficient of viscous damping (lb-sec/ in. ) 

time (sec) 

dimensionless time 

Kt, in free vibration 

oot, in forced vibration 

characteristic root; first five values tabulated in Section II (dimensionless) 

damping factor (dimensionless) 

CiL 


2WK 


, in free vibration 


_ Cg_ 
Wco 


, in forced vibration 
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FREE AND FORCED VIBRATIONS OF CANTILEVER 
BEAMS WITH VISCOUS DAMPING 

SUMMARY 


The fundamental theory for analyzing an undamped vibrating cantilever beam is 
presented in elementary texts on mechanical vibrations. The necessary formulas for a 
complete study of the state of motion and the behavior of a vibrating cantilever beam with 
viscous damping, however, have not been developed and presented for design. 

This repor ^provides the equations of motion for the free and the forced vibrations 
of a cantilever beam with viscous damping. The equations of motion include formulas 
for the bending moment, lateral shearing force, deflection, velocity, and the accelera- 
tion at any desired point of the beam for any chosen time. The coefficient of viscous 
damping is assumed to be constant throughout the length of the beam. ! 

It . is assumed during free vibration that the free end of the beam is initially dis- 
placed some arbitrary distance Y 0 and then released. The mode of vibration and natural 
frequency is defined by the value of a characteristic root. The first five values of this 
root are tabulated) in this report. 

'Jit is assumed during forced vibration that the normally fixed end of the cantilever 
beam is subjected to a lateral displacement of the form y = Y 0 sin cot and that the zero 
slope is maintained. The equations of motion for forced vibration are for the steady-state 
condition only.j 

The equations of motion are presented in dimensionless form in a convenient and 
usable manner. This report also contains the derivations, design curves to evaluate 
some of the parameters, and a set of dimensionless results for each of several sample 
examples. The intent of this report is to bridge the gap between theory and design. 


SECTION I. INTRODUCTION 


Necessary formulas for the state of motion for the free and the forced vibration 
of a cantilever beam are presented in this report. These formulas include necessary 
equations for evaluating bending moment, lateral shearing force, deflection, velocity, 
and acceleration at any desired point along the beam for any desired instant. 



During free vibration it is assumed that the free end of the cantilever beam is 
initially displaced some arbitrary distance Y 0 and then released. Natural frequency 
of the vibrating cantilever beam is dependent upon the mode of vibration, characterized 
by the value of a characteristic root. The first five values of this root are included in 
this report. 

During forced vibration it is assumed that the normally fixed end of the beam is 
subjected to a lateral displacement of the form y = Y Q sin cut, where Y 0 is an arbitrary 
constant and co is the circular frequency of the excitomotor, and the zero slope is main- 
tained. Formulas for forced vibration are for the steady-state conditions only; that is, 
equations which define the state of motion after sufficient time has elapsed for the 
transient terms to decay out of range of significant value. 

Equations of motion are expressed in terms of dimensionless coordinates and 
parameters. These coordinates are axial z and time 0. These parameters are defined 
in the tabulation of the formulas for each type of vibration. One of these parameters is 
the damping factor a which is dependent upon the coefficient of viscous damping. This 
coefficient is assumed constant throughout the length of the beam. Practically any beam 
of uniform cross-section may be analyzed by the procedure presented. 

This procedure and the necessary formulas for each type of vibration are pre- 
sented in a convenient and usable manner in Section II. Design curves for the evaluation 
of the beam parameters and the damping factor are presented in Section III. Also, in 
Section III is a set of dimensionless results for the deflection and bending moment for 
each of several sample cases. The derivation of the formulas is presented in the Ap- 
pendix. It is not necessary to understand this derivation to understand its application. 


SECTION II. APPLICATION 


A. FREE VIBRATION 

1. Procedure . The necessary formulas for the state of motion, bending 
moment, and lateral shearing force in a cantilever beam during free vibration are tabu- 
lated in this section. 

The evaluation of these formulas is straightforward when done in a step-wise 
manner. A brief account of the physical aspects of some of these steps is given in the 
remainder of this section. 

The free vibration of a cantilever beam is characterized by the mode of vibration. 
This mode of vibration is defined by the value of the characteristic root X. The values 
of the first five roots are given in Formula 3. The first two modes of vibration are 
shown in Figure 1 parts (a) and (b) in the definition of symbols. 
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Values of the characteristic roots are solutions of a characteristic equation 
representing four boundary conditions; namely the lateral deflection and slope are zero 
at the fixed end, and the bending moment and lateral shearing force are zero at the free 
end of the beam. 

Only two parameters are required for the initial evaluation of the formulas; the 
damping factor a (Formula 2) and the characteristic root X (Formula 3) . The value of 
the beam parameter K (Formula 1) is used to convert the dimensionless results to 
physical quantities. This parameter ig also included in the definition of the damping 
factor a; . The values of the parameters K and a may be obtained from the design curves 
given in Section in. 

The values of X and a determine the value of the vibration parameter y (Formula 
4) . The significance of this parameter is that the product yK is the natural circular 
frequency of the cantilever beam for the mode of vibration characterized by the value of 
X. The value of y reduces to X 2 when the beam is undamped (a = 0) . 

The integration constant B (Formula 5) is one of two elements of a characteristic 
vector corresponding to the particular value of the characteristic root X. The other ele- 
ment ( integration constant) does not appear in the formulas as it was taken equal to 
one. 

The Z-function (Formula 6) is a function of the dimensionless axial coordinate z 
only. The value of the Z-function at z = 1 is the value of the constant Z t (Formula 7) . 
This constant appears in the T-function (Formula 8) ; a function of dimensionless time 
8 only. The integration constants in the T-function were evaluated in terms of Zj by the 
initial conditions that the displacement is equal to Y 0 and the velocity (dy/dt) is zero at 
the free end of the beam when t = 0. 

The product of the Z- and T-functions is the dimensionless deflection y/Y Q 
(Formula 9) . The remaining results, also dimensionless, are obtained by evaluating 
Formulas 10 through 13. A set of dimensionless results for the deflection and for the 
bending moment (Formulas 9 and 12, respectively) are given in Section III. 

The set of formulas for the free vibration of a cantilever beam are now tabulated 
in the next subsection. The procedure and the formulas for forced vibrations are given 
after this tabulation. 

2. Form ulas . 

Beam parameter 
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Damping factor 


a = 


CfL_ 

2WK 


Characteristic root, X* 

Mode A 

1 1.87510 

2 4.69409 

3 7.85476 

4 10. 99554 

5 14.13717 

Vibration parameter 


( 2 ) 


(3) 


y A 4 - a 2 

Integration constant in Z -function 


(4) 


_ cosh A + cos A 
sinh A + sin A 

The Z -function 

Z = coshA z - cos A z + B(sinh A z - sin A z) 
The Zj constant 

Zj = value of Z -function at z = 1 

= cosh A - cos A + B(sinh A - sin A) 

The T-function 

_ 1 a0 

T = ~ — e (y cos y6 + a sin y6) 

7*1 


(5) 


( 6 ) 


( 7 ) 


( 8 ) 


* A is the root of 1 + cosh A cos A = 0. 
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Deflection, y 



Velocity, V 



- X 4 „ -aO . „ 

— - — Ze sm yQ 
yZt 


Acceleration, A 


A 


-X 4 „ —otQ , 

-yr~ Ze (y cos y6 - oi smyd) 
yZj 


Bending moment, M_ 

M L 2 

— — — = X 2 [ cosh Xz + cos Xz + B(sinhXz + sin Xz) ] T 
EIY q 

Shearing force, Q 


( 9 ) 

( 10 ) 

( 11 ) 

( 12 ) 



sinh Xz - sin Xz + B (cosh Xz + cos Xz)] T 


B. FORCED VIBRATION 


(13) 


1. Procedure . The necessary formulas for the state of motion, bending 
moment, and lateral shearing force during the steady-state condition of the forced vibra- 
tion of a cantilever beam are given’ in this section. The normally fixed end of the canti- 
lever beam is subjected to a lateral displacement y = Y Q sin wt, but the zero slope is 
maintained. A steady-state condition is established when sufficient time has elapsed for 
the transient terms to decay out of range of significant value. 

In the case of forced vibra tion, the beam parameter 0 (Formula 14) is dimen- 
sionless and equivalent to \JuT / K where K is the beam parameter for free vibration 
(Formula 1 in the previous subsection) . The damping factor a for forced vibration 
(Formula 15) is not the same as for free vibration (Formula 2) however, it is still 
dimensionless. 

The parameters (3 and a are the only two parameters required for evaluation of 
a complete set of dimensionless results. The value of the parameter /3 may be obtained 
from the design curves given in Section III. The results are obtained in a step-wise 
manner. The physical aspects of some of these steps and the nomenclature used in this 
presentation are explained in the following discussion. 
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The parameters /3 and a define the four vibration parameters- $, / u, a, and b 
(Formula 16) . The vibration parameters a and b define the symbols si and ti 
(i = 1, 2, 3, 4) which, in turn, define the symbols Sjj and Tjj (Formula 17). The sym- 
bols S^j and Ty are used to evaluate the elements of a four by four matrix (Formula 18) . 

The solution of this matrix equation evaluates the integration constants A, B, C, 
and D and in turn, evaluates the four integration constants E, F, G, and H (Formula 19). 
Various linear combinations of these eight integration constants will be required later. 
The presentation of these combinations is explained below ( in the following paragraphs) . 

The general form of the expression for the dimensionless deflection y/Y Q is as 

follows 


y/Y 0 = Z S sin0 + Z C cos0, 

s c 

where 9 is the dimensionless time cut, Z and Z are functions of the dimensionless 
axial coordinate z only. The significance of the superscripts s and c becomes apparent 
when each is correlated with the sine and the cosine term , respectively. 

g 

The Z function is defined as follows: 

g 

Z = A cosh az cos bz + B sinh az sin bz + 

C sinh az cos bz + D cosh az sin bz + 


E cosh bz cos az + F sinh bz sin az + 

G sinh bz cos az + H sinh bz sin az, 

where the symbols A through H are the eight integration constants evaluated by Formulas 
17 and 18. 

g 

The expression for Z is transformed to the following matrix equation 


s s 

Z = [cosh az sinhaz]M 


cos bz 
sin bz 


+ [cosh bz sinh bz] 1ST 


cos az 


sin az 


s s 

where the eight integration constants are grouped into the matrices M and N as follows 


M S = 

~A 

D~ 

N S = 

~E 

H 


C 

B 


G 

F 


6 


r c s 

The form of the Z function is the same as the form of the Z function, except 

for a different set of eight integration constants. This set is designated by the super- 
script c. The expression for the derivative d n Z/dz n , of any order n, is also the same 
form as the Z-function, except with a different set of eight constants. The second and 
the tliird derivative of Z with respect to z is required for the evaluation of Formulas 25 

and 26. An integer subscript is added to provide the scripted symbols zj, m|, and n|, 

where j is s or c and the integer i designates the order of the derivative with respect to 
z. This scripted notation on the matrices M and N also applies to each of the elements 
within the matrix. 

The non-scripted symbols A through H defined by Formulas 18 and 19 are equiva- 
lent to the corresponding scripted symbols with j = s and i = 0. The value of each of the 
other scripted symbols (that is, the elements of the other matrices M and N) is obtained 
by taking certain linear combinations of the original eight integration constants A through 
H. This linear combination is presented as the linear combination of two, 2 by 2 ma- 
trices (Formula 20). This is illustrated below for the elements Af and G 3 (see matrix 
equation for Mf and N 3 in Formula 20) . 

S s S 

A 2 = A cos20 + B sin2<£ = A^ cos2<£ + B q sin20 , 
and, similarly 

G 3 = F sin3 <£ - E cos3</> = F^ sin3<£ - cos3 <p 

After the evaluation of all the scripted elements (constants) , the desired result 
is obtained by evaluating the appropriate of Formulas 22 through 26; noting that the par- 
ticular z| function is established by substituting the corresponding matrices M| and n| 
into Formula 21. 

The initial results are dimensionless as defined by the expression on the left of 
the equal sign in each of Formulas 22 through 26. The dimensionless time is converted to 
physical time with the relationship 6 = cot. A set of dimensionless values for the deflec- 
tion and bending moment for several examples are presented in Section III. 

The formulas for forced vibration of a cantilever beam are tabulated in the 
remainder of this section. The equations of motion presented in this section are derived 
in the Appendix. 
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2. Formulas. 


Beam parameter, /3 



Damping factor, a 


a - 


CfL 

Wu 


Vibration parameters 


<p = ~tan -1 a 


( 15 ) 


**= 0 ( 1 + a 2 ) 


( 16 ) 

a = p cos $ 
b = p sin <p 

Definition of symbols for the elements of the matrix to evaluate the first four integration 
constants. 


s t = cosh a 

tj = cosh b 

s 2 = sinh a 

t 2 3 sinh b 

s 3 ~ cos b 

t 3 = cos a 

S 4 = sin b 

t 4 = sin a 

Su = s 1 s 3 

Tl3 ~ tjt 3 

S 14 = s lS4 

T 14 = *1*4 

s 23 - s 2 s 3 

T 23 = *2*3 

S 24 = S 2 S 4 

T 24 = *2*4 
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The matrix equation 



Integration constants E, F, G, and H 

E = 1 - A, F = B, G = -D, H = -C. 

Matrices M S , N?, M°, and N C . (i= 0, 2, and 3) 

ill i ' ' 
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M 3 = cos 3 0 


I I + sin 30 I I 

L A D J L B - C J 


G F H -E 

g 

N 3 = -sin 3 0 - cos 30 

EH F -G 


M 3 = cos 30 


D -A 


B -C 


sin 30 


C B 


A D 


( 20 Concluded) 


-HE G F 

N 3 = -sin 30 - cos 30 

-F G EH 


The - function (i = 0, 2, 3; j = s, c)* 

. |"cos bz~"| . ("cos az~j 

z] = [coshaz sinhaz] mH + [cosh bz sinh bz] N? (21) 


sin bz 


sin az 


Deflection, y 


vs c 

dr = Z sin 9 + Z cos 9 
Y o o 

o 


Velocity, V 


Vs c 

— — = Z cos 9 - Z sin 9 
u> Y o o 

o 


Acceleration, A 


£1 s c 

2 " = - (Z sin 9 + Z cos 9) 

or Y ' o o 

o 


* The non-scripted symbols in Formulas 18, 19, and 20 are equivalent to the scripted 
symbols when j - s and i = 0. 


E 


Bending moment, M 

M L 2 s c 

= M (Z 2 sin 0 + Z 2 cos 0) (25) 

o 

Shearing force , Q 

■j§Y" = M 3 (Z® sin 0 + Z° cos 0) (26) 

o 


SECTION IE. SAMPLE RESULTS 


A set of dimensionless values for the deflection and bending moment of a vibrating 
cantilever beam are presented in this section. Design curves for the evaluation of the 
beam parameters K and /3 and for the damping factor a for free vibration are also given. 
The use of these curves is illustrated in each Figure. 

The set of dimensionless results for the cases of free vibration immediately 
follow Figure 4. These cases are for the first two modes of vibration (first two values 
of X) for each of the damping factor a equal to 0 and to 0. 2. These results were obtained 
by evaluating Formulas 9 and 12 of Section II and are presented in Tables 1A through 8A. 

Similar data for the cases of forced vibration are presented after Table 8A. These 
cases are for the beam' parameter (3 equal to 5 and to 10 for each value of the damping 
factor a for forced vibration equal to 0 and 0. 2. These results were obtained by evalu- 
ating Formulas 22 and 25 of Section II and are presented in Tables IB through 8B. 

The format of each table is the same; with the dimensionless time coordinate 0 
(in angular degrees) in the first column and with the dimensionless axial coordinate z 
as columnar headings of the remaining columns. In the case of forced vibration, it is 
recalled that the formulas are for the steady state condition only. That is, the point of 
zero time (0 = 0) is merely a beginning reference point of a repetitious cycle. This 
explains why the results are not zero everywhere along the beam at 0 = 0 when the beam 
is subjected to viscous damping (cv ± 0) . The derivation of the formulas in Section II is 
given in the Appendix which follows Table 8B. 


il 



7 


Data, for illustrative example: 
El = 1. 6 x 10* lb-in. 2 
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FIGURE 2. THE BEAM PARAMETER K 





FIGURE 3. THE BEAM PARAMETER j3 
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TABLE 1A. DEFLECTION OF CANTILEVER BEAM - FREE VIBRATION 
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TABLE 2B. DEFLECTION OF CANTILEVER BEAM - FORCED VIBRATION 
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330 -.500 -.512 -.510 -.463 -.368 -.245 
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TABLE 3B. DEFLECTION OF CANTILEVER BEAM - FORCED VIBRATION 
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TABLE 4B. DEFLECTION OF CANTILEVER BEAM - FORCED VIBRATION 
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TABLE 5B. BENDING MOMENT IN CANTILEVER BEAM - FORCED VIBRATION 
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TABLE 8B. BENDING MOMENT IN CANTILEVER BEAM - FORCED VIBRATION 
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APPENDIX. DERIVATIONS 


A. INTRODUCTION 


The derivations of the formulas presented in Section II are given in this sec- 
tion. These derivations stem from the solution of a differential equation whose develop- 
ment for undamped vibration is presented in texts on elementary strength of materials 
and mechanical vibrations. The presentation in this report treats primarily the develop- 
ment of the formulas from the solution of this differential equation. 


B. INITIAL DIFFERENTIAL EQUATION 

A brief account of the development of the initial differential equation for a 
vibrating cantilever beam is given in this section. This development begins with the 
equation of equilibrium and the moment-curvature relationship written below. * 


1. Equilibrium Equations 


dx 


m 


if. . 

at* 


at 


™ = Q (i) 

ax 


2. Moment-Curvature Relationship 


M 


* E I 


9?y 

8x 2 


( 2 ) 


In the above relationships, Q is the lateral shearing force and M_is the bending 
moment in the beam. The symbols m and c are the mass and the viscous damping coef- 
ficient per unit length of the beam, respectively. 


The initial differential equation is developed by eliminating M_ and Q from Equa- 
tions 1 and 2. This yields 


E I 



*2- _ C *L 
at? at 


(3) 


The formulas for each type of vibration are developed from a solution of this dif- 
ferential equation ( Eq. 3) . The first step in this development is the transformation of 
Equation 3 to dimensionless form. This transformation is slightly different for each type 
of vibration. This transformation and the subsequent derivations for each type of vibra- 
tion is given in separate subsections. The formulas for free vibration are developed in 
the next subsection. 


* These relationships for the undamped vibrations of a cantilever beam are developed 
and presented by W. T. Thomson, "Mechanical Vibrations," Prentice-Hall, Inc. 1948, 
Chapter 6. In this report, the positive direction of the coordinate y is changed. 
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C. FREE VIBRATION 


The derivation of the formulas for free vibration is given in this section. 
This derivation begins by introducing the following coordinate transformations and 
symbol definitions. 


x = zL, 



c 


C 

L 


K 2 = 


Elg 

WL 3 


Q! = 


2WK 


( 4 ) 


e = Kt. 


The above definition of K and a is the same as Formulas i 


The definitions and coordinate transformations in Equation 
3 to the following form. 


<*L + <dL +2a <£L = 

q-A q/i 2 " ao 


80 


de 


o. 


and 2 of Section II. 

4 transform Equation 

(5) 


The solution of Equation 5 is obtained by the method of separation of variables, 
that is, the deflection y is taken to be of the form 

y/Y Q - ZT, (6) 


where Z is a function of the dimensionless axial coordinate z only, T is a function of 
the dimensionless time 6 only and'Y o is an arbitrary constant. 


Substituting Equation 6 into Equation 5 and separating the variables yields the 
following two ordinary differential .equations. 


£-x*z-o. 

dz* 


4?? + 2 a ~ + A 4 T - 0, 


dr 


de 


(a) 


(b) 


(7) 


where X is a constant, which will be evaluated later from the specific conditions imposed 
upon the vibrating beam . 


Equation 7a is written below in a form which satisfies the boundary conditions 
that the deflection and slope are zero at the fixed end of the cantilever beam; in other 
words , z = dZ/dz = 0 at z = 0. 


Z = A(cosh X z - cos X z) + B(sinh X z - sin X z). 


( 8 ) 
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The remaining two boundary conditions are that the bending moment and lateral ’ 
shearing force are zero at the free end of the beam. This corresponds to setting the 
second and the third derivatives of Z with respect to z equal to zero at z = i. Thege 
conditions yield the following matrix equation. 


— 

cosh A + cos A 

sinh A + sin A 


A" 

sinh A - sin A 

cosh + cos A 


B 


(9) 


A non-trival solution of Equation 9 is obtained only when A is the characteristic 
root of the coefficient matrix. Equating the deteriment of this coefficient matrix to zero, 
yields the following expression for the evaluation of the characteristic root. 

1 + cosh A cos A = 0. (10) 

The first five roots of Equation 10 are tabulated in Section II as Formula 3. 
Accuracy greater than eight significant digits is required to obtain the higher roots be- 
cause of the characteristic behavior of Equation 10 at these higher values. The mode of 
vibration is characterized by the value of the characteristic root, A. The first two modes 
are shown in Figure la and lb in the definition of symbols. 

For each characteristic root there corresponds a characteristic vector whose 
elements are the integration constants A and B. The value of the integration constant A 
is arbitrarily taken equal to one. The integration constant B is evaluated from the first 
of Equation 9; this yields: 


„ cosh A + cos A ,, 

B = - — ~ - : — . (11) 

sinh A + sin A 

Equation 11 is Formula 5 of Section n. The evaluation of A and B completes the 
definition of the Z-function in Equation 8 and is Formula 6 of Section II. The constant 
Zj (Formula 7 of Section II) is the value of the Z-function at z = 1 and will be used later. 
This discussion continues with the solution of Equation 7b. 

The general solution of Equation 7b is 

q/Q 

T = e~ (E cos y 0 + F sin y 0) , ( 12) 

where E and F are integration constants, and y =f\j A 4 - a 2 . During this solution it is 
assumed that A 4 is greater than a 2 which must be true in order to have oscillatory motion. 

It is assumed that the initial conditions for free vibration are that the deflection 
is equal to an arbitrary displacement Y Q and that the velocity is equal to zero at the free 
end of the beam when t — 0; mathematically, 
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y =y q . 


when 0 = 0 and z = 1. 


( 13 ) 


and 


djL = 
d0 


0, 


Imposing these conditions upon Equation 6 yields 
E = 1/Zt, and F = a/y Z t . 


(14) 


Substituting Equations 14 into Equation 12 completes the definition of the T-f unction; 

that is, 


T = 


-aQ 


(y cosy 0 + a sin y0) , 


which is Formula 8 of Section II. 


(15) 


Introducing the above definitions of the Z-and the T-functions into Equation 6 com- 
pletes the derivation of Formula 9 of Section II. The derivation of the remaining formulas 
for the free vibration of a cantilever beam (Formulas 10 through 13) are developed in a 
straightforward manner by taking the appropriate derivative of Formula 9. The defini- 
tions and the corresponding coordinate transformations for this are tabulated below. 


Velocity 


V -9£_ k 93L_ ky ^ — 
V ~ at ~ K 90 ~ K Y 0 z d0 


Acceleration 




_ 

at 2 


90* " o d0 

Bending Moment 

_ F T EI Y o d!z 


M 


ax 2 l 2 az 2 


dz* 


(16) 


Shearing Force 

Q = EIt^- 


8x J 


_ E I 8 J y 

l 3 az 3 


EI Y o d*Z 


dz 


Formulas 10 through 13 are verified by substituting and performing the operations 
indicated in Equations 16. These operations are elementary and straightforward. 
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This concludes the discussion on the derivation of the formulas for the free vibra- 
tion of a cantilever beam presented in Section II. The formulas for the forced vibration 
of a cantilever beam is given in the next subsection. 

D. FORCED VIBRATION 


The formulas for the forced vibration of a cantilever beam are derived in this 
section. This derivation begins with the initial differential Equation 3. 

In the case of forced vibration, it is assumed that the normally fixed end of the 
cantilever beam is subjected to an exciting displacement of the form y = Y Q sin cot, where 
Y 0 is a constant and co is the circular frequency of the excitomotor. The formulas pre- 
sented in Section II for this type of vibration pertain only for the steady state condition. 
This condition is characterized mathematically by a particular rather than by the general 
solution of Equation 3. 


The initial step for the derivation of this particular solution is the transformation 
of Equation 3 to dimensionless form, by introducing the following symbols and coordinate 
transformations . 


x = zL, 



0 = cot, 



A _ WL 3 CO 2 
P Elg K 2 ’ 


a 


CK 

Wco 


(17) 


The above definitions for /3 and a are Formulas 14 and 15 of Section II, respec- 
tively. 


Imposing these definitions and transformations upon Equation 3 yields 


9 V 

dz 


+ /J 4 (°5t- +a SsS 
p ya0 2 de 


= o. 


(18) 


A particular solution of Equation 18 is assumed to be of the form 

y/Y Q = Z S sin 0 + Z° cos 0, 
s c 

where Z and Z are functions of the axial coordinate z only and Y Q is a constant. 


(19) 


Equation 19 is a solution of Equation 18 provided that the Z -functions satisfy the 
following two fourth-order ordinary differential equations simultaneously. 
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( 20 ) 


~ ? (Z S + a Z°) =0, 

^ (Z c - 0!Z S ) = 0. 

The following eighth-order differential equation is developed by combining Equa- 
tions 20 to eliminate Z c . 


d 8 Z S 

dz 8 


d 4 Z S 


“ 2 Z 3 * dz 3 — + ^( i + a2 ) z = °- 


( 21 ) 


21 . 


The following symbols are introduced for writing the general solution of Equation 


/? ( 


1 + a 




a 


a = n cos <p 


b = n sin c p. 


( 22 ) 


The definitions in Equations 22 are Formulas 16 and are determined by the roots 
of the auxiliary equation of the eigntn-order differential Equation 21. The general solu- 
tion of Equation 21, in matrix form is 


Z = [cosh az sinh az] M 1- 


where the eight integration constants 
matrices M s and N s . 


— — 

cos bz 

s 

+ [ cosh bz sinh bz] N 

cos az 




(23) 

sin bz 


sin az 



are presented as the elements of the following two 


M S = 


A D 
C B 


N = 


E H 
G F 


(24) 


where the single superscript on the right is applicable to each of the elements in the 
matrix. 


The form of the Z° function is identical to the form in Equation 23, except for a 
different set of eight integration constants. This set is designated with the superscript 
c. In order to satisfy Equations 20 simultaneously, one set of the eight integration con- 
stants must be dependent upon the other set. To show this and then to evaluate the re- 
maining eight integration constants it will become necessary to develop the first four 
derivatives of the Z -functions. 
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The form of the derivative of the Z-function, of any order, is the same form as 
Equation 23, except with a different set of eight constants. These other sets of constants 
and the order of the derivative are designated with an integer subscript on the matrices 

M and N. This creates the scripted symbols zj, m|, and n|, where j is s or c and 

i = 0, 1, 2, 3, or 4. The non-scripted symbols developed above are equivalent to the 
scripted symbols when j - s and i = 0. The superscripts s and c are omitted when the 
relationship is for both as is the case in the following discussion. 

The elements of the sequential matrices (corresponding to the sequential deriva- 
tives) are obtained from a recurrence formula. Two transformation matrices, desig- 
nated with the symbols J and K, are introduced for writing this recurrence formula. 



0 1 


r- — i 

0 -1 

J = 


, and K = 



_1 0 


1 0 


The recurrence formulas for the sequential matrices M and N are written below. 


M. = aJM. 4 + bM. K 

i i-i i-i 

(i = 1, 2, 3, 4) 

N. = bJN. + aN. K 

l i-i i-i 


(26) 


The elements of the matrices for i ^ 2 are expressed in terms of the matrices M Q 
and N Q by a sequential substitution process. This process is illustrated for the first of 
Equation 26 with i = 1 and 2. 

M< = aJM + bM K , 
o o 

(27) 

M 2 = aJMi + bM t K . 

Before substituting the first into the second of Equation 27, the following proper- 
tv of matrices J and K is recognized. 

JJ = - KK = I, (28) 

where I is the identity matrix of order two. 

Accordingly, 

M 2 - (a 2 - b 2 ) M + 2abJM K, 

which, upon substituting for a and b from Equations 22 and using known trigonometric 
identities, transforms to the following form 
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(29) 


M 2 = i u 2 (cos 2<p + sin 2 <f> JM^K) . 

The remaining matrices are developed in a similar manner. The final form of 
matrices M i and Nj (i = 1, 2, 3, 4) are tabulated below. 

M * = p (cos cf> JM q + sin (p M Q K) 

M 2 = p 2 (cos 2<p + sin 2<p JM q K) 

M 3 = p 3 ( cos 3(f> JM q + sin 3 <f> M q K) 

M 4 = p 4 (cos 4(f) M q + sin 4 (p JM q K) 

(30) 

Nj = jli (sin <p JN q + cos (p N^K) 

N 2 = p 2 (-cos 2(f> N q + sin 2 <f> JN^K) 

N 3 = -p 3 (sin 3 (p JN^ + cos 3 cp N q K) 

N 4 = p 4 (cos 4(p N q - sin 4<p JN q K) . 

Any derivative of the Z-function is obtained by substituting the appropriate matrix 
M and N of Equations 30 into 23. The order of the derivative is designated by the integer 
subscript. 

The above relationships are used to evaluate the two sets of eight integration con- 
stants; that is, the s-set and the c-set. The eight integration constants in each set are 
designated by the capital letters A through H and are presented as the elements of the 
matrices M and N; each with its corresponding superscript. The location of these con- 
stants within the matrices is as shown in Equations 24. 

The integration constants in the c-set are expressed in terms of the s-set by 
substituting Equation 23, with the appropriate definition in Equation 30, into the first of 
Equations 20. This yields 

M C = JM S K , and N C - -JN S K, (31) 

00 00 ' 

or, in expanded form 
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The set of matrix equations in Formulas 20 in Section II is verified by substituting 
Equation 31 into the appropriate definition in Equations 30 and expanding the matrices 
and upon noting that the p 2 and p 3 terms are omitted in Formulas 20 as they are included 
in Formulas 25 and 26. 

The eight integration constants in the set are evaluated by applying a set of four 
boundary conditions to each of the two Z-functions. These conditions are tabulated be- 
low. 



at z = 0 


at z = 1 


(a) 


(32) 


(b) 


The first two of Equations 32a ensure that the displacement at z = 0 is of the 
form y = Y q sin0, where Y 0 is the amplitude of the exciting displacement, (see Equation 
19) . The second two of Equations 32a ensure that the slope dy/dx is zero at z = 0. 
Equations 32b ensure that the bending moment and lateral shearing force are zero at 
z = 1. 


The value of the Z-function at z = 0 (Eq. 23) reduces to 
Z = A + E. 


Equations 32a transform to the following set of equations. 


A S + E S = 1 
o o 

A c _c 
A + E = 0 
o o 

A® + E® = 0 


A° + E° = 0 


(33) 


( 34 ) 
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Equation 34 are expressed in terms of the s-set of integration constants by 
using Equation 31 and the appropriate definitions in Equations 24, 25, and 30. This 
result is written below. 

“\ 


A S + E S = 1, 
o o 


B S - F S = 0, 
o o 


(a) 


cos 0 (C® + H®) + sin0 (D® + G^) 
-siiicp (C® + H®) + cos0 (D® + G®) 
or, alternatively for Equations 35b. 



(b) m 


C s + H s 
o o 


D S + G S = 0 . 
o o 


(C) 


The following symbols are introduced for writing the value of the Z -function at 
z = 1 , and are used in developing the expression for the remaining four boundary condi- 
tions in Equations 32b. 


Sj = cosh a 

t 4 = cosh b 

s 2 = sinh a 

t 2 = sinh b 

s 3 = cos b 

t 3 = cos a 

s 4 = sin b 

t 4 = sin a 

S 13 = s 4 s 3 

T 13 = *1*3 

S 14 - Sl s 4 

T 14 = *1*4 

®23 = S 2 S 3 

T 23 = t 2 t 3 

S 24 = S 2 S 4 

T 24 = t 2 t 4 


The definitions in Equations 36 are the same as those of Formulas 17 in Section 

II. 
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The value of the Z-function at z = 1 is written below in terms of -the above sym- 
bols; first in matrix form and then in algebraic form. 



s 3 “ 



— — 

^3 

z = [s 4 s^ M 


+ [ tj t 2 ] 

N 



s 4 





Z = AS 13 + BS 24 + CS 23 + DS 14 + ET 13 + FT 24 + GT 23 + HT 14 


Equations 37 are valid for any set of the scripted symbols introduced previously. 
In other words, the expression for any Z-function and any of its derivative is obtained by 
merely applying a consistent set of scripts to each of the symbols Z and A through H in 
Equations 37. To illustrate, the second derivative of Z s and Z c are written below. 


z = [ Sl s 2 ] m; 


“s 3 " 

c 

t 3 


+ [ti k] N 


_s 4 


_t 4 _ 


+ [ti t 2 ] Ng 

where, upon substituting Equation 31 into the 
then simplifying with Equation 28, 

= p 2 (cos + sin 2<pJ M^K) 

Mg = p 2 (-sin 2<£M^ + cos 2<p J M^K) 

Ng = p 2 (-cos 2<£N^ + sin 2<pJ N^K) 

N 2 = p 2 (sin 2(p + cos 2<p J N^K) . 


z 2 ~ [ s i s 2l M g 


s 3 


s 4 


^3 

U 


(38) 


appropriate definition of Equation 30 and 


(39) 


The algebraic representation of the first two conditions in Equations 32b are ob- 
tained by substituting the appropriate set in Equations 39 into Equations 38. These initial 
algebraic equations are simplified by a large degree by taking certain linear combinations. 
These linear combinations are illustrated below for Equations 38. 


sc s 

cos 2 <p Z 2 - sin 2<t> Z 2 = [s 4 s 2 ] M q 


r -i 

S 3 


- [t 4 tj N 


sin 2 <p zf + cos 2<p Z 2 - [ s 4 s 2 ] J M^K 


^3 


s 3 

L -J 

s 

~h 

S 4 

+ [ti t 2 ] jn q k 

h 


( 40 ) 
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The corresponding relationships for the last two of Equations 32b are developed 
in a similar manner and are written below. 


SC s 

cos 30 Z 3 - sin 30 Z 3 = [sj s 2 ] J M q 


SC s 

sin 30 Z 3 + cos 30 Z 3 = [sj s 2 ] M q K 


®3 

s 4 

s 3 

S 4 


- [t t t 2 ] N K 
o 


- ttj t 2 ] J N*" 


*3 

t 4 

*t 3 


1>J 


(41) 


Equations 40 and 41 are transformed to the matrix equation written below by sub- 
stituting Equations 24 for matrices and N® and introducing the symbols Sy and as 
defined in Equations 36. 


S 13 

S 24 

S 23 

Sl4 


A S 

0 


CO 

H 

1 

T24 

T 23 

t 14 


V 

0 


0 

_S 24 

S 13 

-S i4 

S 23 


B S 

o 


T24 

-T 13 

Tm 

-T 2 3 


F S 

0 


0 

S 23 

Sj 4 

s i3 

S 24 


c s 

0 


-Tl4 

T 23 

-T24 

Tl3 


G S 

0 


0 

1 

1 

CO 

£ 

S 23 

-S 24 

S 13 


D S 

0 


T 23 

Tl4 

T i3 

T 24 


H S 

o 


0 


(42) 


The elements of the second vector are related to the elements of the first vector 
in Equation 42 by the following expressions (see Equations 35a and 35c). 


E = 1-A, F = B, G = -D, and H = -C. 


(43) 


The scripted notation is no longer required. Substituting Equations 43 for the 
elements of the second vector, transforms Equation 42 to the following four by four 
matrix. 




— 


r— — 

®13 + T 13 s 24 - HiA, S 23 + T 14 S 14 + T 23 


A 


T» 

-S 24 + T m S 13 + T 13 -S 14 - T 23 Sg 3 + t 14 


B 


T 24 

S 23 “ T 14 S 14 “ T 23 S 13 + T 13 S 24 “ T 24 


C 


"Ti4 

“ S i4 + T 2 3 S 23 - T 14 -S 24 + Ta S 13 + T 13 


D 


T 23 




_ _ 


■— — 


(44) 


Equations 43 and 44 are Formulas 19 and 18 in Section II, respectively. This 
verifies each of the Formulas 14 through 22 of Section II. The remaining formulas (23 
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through 26) are verified by taking appropriate derivatives of the expression for the beam 
deflection (Eq. 19) . This is done in a straightforward manner by substituting and per- 
forming the operations indicated in the following definitions. 

Velocity 


V = - 7 ^- = co 7 ^— = co Y (Z S cos 6 - Z C sin 0 ) , 
8t 90 o v o o ' 


Acceleration 


A = = a ) 2 - - u ) 2 Y (Z S sin 0 + Z° cos 0) , 

— dr 86 n « ' 


o o 


Bending Moment 


d\ El EI Y 0 S C 

M = E I ^2 = ~jj~ ( Z 2 sin 6 + Z 2 cos 0 ) • 


Shearing Force 


, , EI Y 

0 3 v EI 8 3 v o „s c 

Q = El ~^ = I? d^^~lT (Z 3 Sin0 + Z 3 cos *> 


(45) 


The first two of Equations 45 are Formulas 23 and 24, respectively. Formulas 
25 and 26 are verified with the last two of Equations 45 by recognizing that the p n terms 
appearing in Equations 30 are omitted in Formulas 20 but are included in Formulas 25 
and 26. 


George C. Marshall Space Flight Center 

National Aeronautics and Space Administration 
Huntsville, Alabama, May 13, 1965 
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